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^5 ■ Abstract 
< 

I We study extensions of the Standard Model where the charged current weak 

interactions are governed by the CKM matrix and where ah tree-level decays 
^ ' are dominated by their Standard Model contribution. We constrain both 

o\ ' 

fS) ' New Physics and the Standard Model contributions to the mixing amplitude 

m 
o 



of the neutral B system using the experimental results on Ru-, Amd,s, 
and a^Ks- We present new results concerning models with minimal flavor 
violation and update the relevant parameter space. We also study the left- 
Q_i' right symmetric model with spontaneously broken CP, probing the viability 

of this model in view of the recent results for a^Ks ^'^cl other observables. 
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\^ ■ I. INTRODUCTION 

For more than 30 years the only unambiguous indication for CP violation (CPV) was 
found in the neutral kaon system. Even though the CP violating parameter 

Ek = (2.280 ± 0.013) X 10"^e''^^ with 0^ ~ 7r/4 (1) 

has been measured rather accurately hadronic uncertainties obscure the determination 
of the fundamental CKM parameters 0. However, there is growing evidence that CPV 
also occurs in the B meson system. Specifically, indications for a CP asymmetry in Bd 
J/^Ks decays have been reported by the ALEPH §, BaBar |]], BELLE §, CDF § and 
OPAL collaborations. Combining all the measurement gives: 

a^^^ = 0.51 ±0.18. (2) 
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It is essential to understand the impact of the present and upcoming experimental data 
on the SM and its extensions. In this paper we focus on a large class of NP models where 
the charged current weak interactions are governed by a 3 x 3 unitary matrix, just like 
in the CKM picture. Moreover, we make the reasonable assumption that all tree-level 
decays are dominated by their SM contribution. However, flavor changing neutral current 
processes are sensitive to New Physics (NP) because, in the SM, they appear only at the 
loop level. 

After a brief update on the SM picture, we study in detail, both analytically and 
numerically, the allowed parameter space of the ratio and the phase difference between the 
NP and the SM contributions to the mixing amplitude of the neutral B system. We use the 
experimental results on Arud^s, Ek and a^Xg fo'^ the generic NP models. Besides the 
two general assumptions mentioned above we introduce additional constraints that apply 
to specific NP models. First, we assume that the ratio between the mass differences of the 
Bfi and the Eg neural meson system is unaffected by NP, which is valid, e.g., in models 
with Minimal Flavor Violation (MFV). Second, we impose a small CP violating phase, as 
for example in left-right symmetric (LRS) models with spontaneous CPV. 

We devote special attention to models of MFV. We show that the three parameters p, 
7] and Ftt can be obtained analytically from the experimental results on Am^ and Ek- 
Additional constraints on the parameter space arise from the measurement of a.^Ks and 
the upper bound on Arria/ Amg. 

Finally, we investigate the impact of the uncertainty of various parameters in LRS 
models with spontaneous CPV on recent claims that these models would be ruled out by 
a significant CP asymmetry in Bd ^ J/ip Kg decays. 



II. STANDARD MODEL PICTURE 



CPV appears naturally in the Standard Model with three generations of quarks. It 
can be attributed to a single phase, 5km i appearing in the CKM matrix that describes 
the weak charged- current interactions. Among the four Wolfenstein parameters (A, A, p, if) 
that provide a convenient description of the CKM matrix only two have been determined 
with a good accuracy [Q, namely 

A = 0.2237 ±0.0033 and A = 0.819 ± 0.040 . (3) 

The other two parameters (p, r/) describing the apex of the unitarity triangle, that follows 
from the unitarity relation VudV*fj + VcdVch + '^tdYth — i have a rather large uncertainty. As 
a consequence the CP violating phase 5km = arctan(r7/p) of the SM is relatively weakly 
constrained. 

In the following we review briefly how the allowed region for the two Wolfenstein pa- 
rameters p and ri is obtained (for more details see e.g. [0,|Tn| and references therein). At 
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present these parameters are experimentally constrained by the measurement of 
the observed — mixing parameter Am^, the lower bound on the B^ — B^ mixing 
parameter Am<(, the measurement of Ek and the combined result for a^Ks- Note that for 
the purpose of this paper it is sufficient to work to leading order in and we therefore 
neglect all subdominant contributions. 

The distance from the origin to the apex (p, rf) of the rescaled unitarity triangle, 



(4) 



is proportional to iKib/Kfel, which dominates its uncertainty. The ratio can be 

determined from various tree-level B decays (see e.g. 0] and Tab 1.), restricting Ru to 
the interval 



Ru e [0.34, 0.43] 



(5) 



The distance between the apex of the rescaled unitarity triangle and the point (1, 0) in 
the p — f] plane is given by: 



VtlVu 



(6) 



Rt can be determined from the experimental results for mass differences Arriq in the B^~B^ 
(g = d, s) systems by means of the SM predictions (see e.g. [p|,p!0|): 



q2 

^VBrnB^m^wflBB^Soixt) ■ \VtqVt 



*|2 
th\ 



(7) 



Here Gp is the Fermi constant, 77^ is a QCD correction factor, m^^ is the 5° mass, my^ 
is the ly-boson mass, /^^ is the Bq decay constant, Bb^ parameterizes the value of the 
relevant hadronic matrix element and 5*0 (xj) gives the electroweak loop contribution as a 
function of the top quark mass = Xtvri^. Also the ratio Amd/ lS.ms is very useful since 
it is simply related to Rt via 



Am^ 
Am, 



p2^ ^B, 

72 

r ^B, 



(8) 



where ^ = /bJ fs^ ■ \J BbJ Bb^ ■ 

The measurement of the Ek parameter in combination with the SM prediction deter- 
mines a region between two hyperbolae in the p — rj plane, described by (see e.g. P,p!0|): 



V 



;i - p)v2So{xt)A^ + P, A^Bk ^ 0.226 . 



(9) 



Here 772 is a QCD correction factor and Pc characterizes the contributions with one or two 
intermediate charm-quarks. 
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Finally the recent results for the CP asymmetry in the ^ J/ip Ks decay in eq. (0) 
imply a further constraint in the p — r] plane via the SM prediction 

a^^Ks = sin 2/3, (10) 

where 

In terms of the Wolfenstein parameters we have 

Combining the above formulae and using the input parameters in Tab. 1 we obtain the 
allowed region for (p, r]) shown in Fig. |I|. Note that we naively combine the aforementioned 
constraints in order to determine the allowed region. (For more sophisticated analyses see 
e.g. refs. Il|JT|PJT^ 



III. GENERIC NEW PHYSICS MODELS WITH UNITARY MIXING MATRIX 

We focus our analysis on a large class of new physics (NP) models with the following 
two features (c.f. refs |T5|-[T7| and references therein): 

(i) The 3x3 mixing matrix for the charged current weak interactions is unitary. 

(ii) Tree-level decays are dominated by the SM contributions. 

Note that assumption (i) holds in all SM extensions with three quark generations. As- 
sumption (ii) implies in particular that the phase of the B — ipKs decay amplitude is 
given by the Standard Model CKM phase, aYg{VcbV*g) . Both assumptions are satisfied by 
a wide class of models, like supersymmetry with i?-parity, various LRS models, models with 
minimal flavor violation (MFV) and several multi-Higgs models. Since some of the models 
mentioned above have additional features we investigate the following three scenarios: 

(a) Only assumptions (i) and (ii) are valid. 

(b) Assumptions (i) and (ii) apply and Arrid/ Arris is unaffected by the presence of NP. 

(c) Assumptions (i) and (ii) apply and 6km is small. 
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Scenario (b) contains, for example, certain multi-Higgs models with "flavor conserva- 
tion" and models with MFV [p!0|,p!8|-p0| . Scenario (c) applies to models of approximate 



CP symmetry and LRS models with spontaneous CPV [pTHzq]. To be specific we take 
I ta,n{6KM)\ < tan(5™'"^) = 0.25 for these scenarios. 

The NP effects relevant to our analysis modify the Bd—Bd mixing amplitude, M12 — |ri2. 
Assumption (ii) implies that the impact of NP on the absorptive part of this amplitude is 
negligible: 



\2 ^ Tff . (13) 



The modification of M12 can be parameterized as follows (see refs. |T^|T^J2^ and refs. 
therein) : 

-rle'^'^. (14) 



The variables and 6d relate directly to the experimental observables. The results for 
Anid^s constrain r^, while the recent measurement of a^Xs restricts the value of the phase 
26d. However, to investigate a specific model it is convenient to use different variables, h and 



0", and describe the NP mo difications to the amplitude M12 = Mg^ + Mf/ via 



The variables h and a are related to and 6d through: 

rje^"^'' = l + he''' . (16) 

Our goal is to determine the allowed parameter space for h and a. To this end we first 
constrain and 29d using the relevant experimental results on R^, Arrid^s, and a^Ks- 
Then we show how to translate the ranges for and 26d into the allowed region in the 
h — a plane. All the relevant parameters we use in our analysis are summarized in Tab. 1. 
For most of these parameters we adopt the values of ref. 0]. For simplicity we scan over 
the ranges indicated in Tab. 1 both for the theoretical and experimental parameters. 



A. Constraining and 6d 



The parameter is constrained by the results for Amd,s [0)0 • models that satisfy 
the assumptions (i) and (ii) the expressions for Arud^s in eq. (|^) have to be modified as 
follows 



q2 

= -^VBmB,ml,fl^BB,Soixt) \VtgV;f ■ rl , (17) 
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where q = d,s and is the analog to for the Bg system. Moreover Am^^/Am^^ satisfies: 



'■ ■ H 
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which reduces to the SM expression in eq. when = r^. Note, that the constraint on 
Ru given in eq. (|) is unchanged by NP, since it is dominated by tree-level amplitudes. 
Using the values in Tab. 1, from eq. (|T^ it follows that: 

0.69 <rd- Rt{p,v) < 1-03, (19) 

while eq. (0) implies that 

-Rt{p,r]) < 0.99. (20) 



The allowed range for the phase 9d is restricted by the measurement of a^Xs IHI- Due 
to the two-fold ambiguity of the sin function, the one-sigma result for a^Ks i^i eq. (^ 
corresponds to two allowed intervals for 2{9d + 13) ||T5| , [T6| : 



19° < 2(0d + /3) < 44° or 19° < 180° - 2(0^ + /3) < 44° . (21) 

We use the above results for two types of analysis: In our numerical analysis we scan over 
the p — rj plane and calculate the allowed intervals for according to eq. (|19D [and eq. ( ^Of ) 
in scenario (b)] and 9d from eq. (^) for all points (p, 77) that fulfill the constraints relevant 
to a particular scenario (c.f. Fig. |^ for scenario (a). Fig. ^ for scenario (b) and Fig. ^ for 
scenario (c), respectively). 

To get better insight into the numerical results we have performed an analytical calcu- 
lation, where we use the maximal allowed region for and 9d that follow from the global 
constraints on Rt and [3. Such an approach is particularly useful, because the analytical 
constraints can easily be updated using the forthcoming experimental data. Moreover the 
system of constraints may be extended for different classes of NP models. For scenario (a) 



we get (see appendix |A 1| for details): 



Td E [0.5, 1.8] and 29d E [-32°, 212°] . (22) 

Note that the interval for 29d in eq. (^) is independent of the upper bound on a^Xs long 
as {29d)maxi > (2^^rf)™"2. With the upper bound on Ru in eq. @) two distinct intervals for 
29d only appear if {a^Ks)max < 0.63. This explains why the global constraint on 29^ we 
obtain from the measurement of a^Xs eq. (0) almost coincides with the one in ref . ||TB[ , 
even though our upper bound on a^Xs is much more stringent. 

For scenario (b) the lower bound on is somewhat stronger, since from eq. ( pOf ) it 
follows that Rt{p, Tj) < 0.99, which gives > 0.7. In scenario (c) we have 
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rd e [0.5, 0.8] with 29^ G [12°, 51°] U [129°, 168°] , 



(23) 



for the allowed regions with p < and 

ra G [1.0, 1.8] with 2^^ G [-1°, 64°] U [116°, 181°] , (24) 
for the allowed regions with p > (c.f. Fig. |). 

B. Constraining h and a 

Our first observation is that from the maximal value for one obtains an upper bound 
on h. Eq. (|T6|) implies that 



h{rd, Od) = ^rj-2rjcos2ed + l < (rr")' + 1 ~ 4.24 , (25) 

where the upper bound arises from setting cos 26 d = — 1 and = r^°^. Still, one can 
obtain much better bounds using the correlation between h and a. Here we only sketch 
how to translate the constraints one obtains for and 26d into those for h and a. The 
details can be found in appendix [A 2\ The basic idea is simple. We solve eq. (|16D for h, 
once eliminating 6d, 



h^{(J, Td) = — cos cr ± \J cos^ a + rj — 1 , (26) 
and once eliminating r^, 

, ±. . . sin2grf sm{29d±(T) , . 

h fT, Od) = — — r^— — . 27 

sm a — sm 26d 

Using the above equations the ranges for and 6d can be translated into a region in the 
h — a plane. The overlap of the two areas from eq. (|2^) and eq. ( P?] ) corresponds to the 
allowed region in the h — a parameter space. While this method is straightforward one has 
to treat carefully the discrete ambiguities of the trigonometric functions (see appendix). 

Using eq. ( ^B]) and eq. ( P7| ) the global constraints on and 9d in eqs. (P^)-(PID can be 
readily translated into the allowed region for h and a. Combining all the resulting intervals 
one obtains the light areas, whose boundaries are given by h^{a,rd) and h^{a,9d) with 
and 9d at the edges of their allowed intervals (see appendix [A2| for more details). 

In order to obtain the exact allowed region in the h — a plane we have performed a 
numerical analysis proceeding as follows: Rather than using the global constraints on 
and 9d corresponding to the entire allowed region in the p—r) plane, we compute the possible 
ranges for and 9d for each permissible value of p and rj. That is we scan over the p — rj 
plane and for all points (p, r]) that fulfill the constraints relevant to a particular scenario, 
we calculate the allowed intervals for according to eq. (|T9|) [and eq. ( [20| ) in scenario (b)] 



7 



and 9d from eq. (0). The resulting intervals are then translated into the corresponding 
region in the h — a plane, just as we did in the "global" analysis, but for each point (p, 77) 
separately. Combining all the resulting intervals one obtains the allowed regions indicated 
as the dark areas in Fig. I, Fig. I and Fig. 0. They are contained within the light areas 
resulting from the analytic boundaries from the global constraints. Even though in general 
the numerical results are not much more constraining than the global ones, there exist 
some regions where the numerical constraints are significantly stronger than the analytical 
boundaries. For example, there appear "holes" in the dark gray regions of Fig. ^ and Fig. |^ 
that correspond to the excluded region within the Ru annulus, which has been ignored in 
"global" analysis (light gray regions). Even more important, the lower bound on h in 
scenario (c) (see Fig. ^ from the numerical analysis turns out to be significantly stronger 
than in the global analysis, which is marginally consistent with /i = 0. 



C. Summary of Results 



(a) Generic models: All models that satisfy the assumptions (i) and (ii) must respect 
the condition in eq. (^), which corresponds to an annulus around the origin in the 
p — 1] plane as shown in Fig. ^. The resulting allowed region in the h — a plane is 
presented in Fig. |^. The admissible range for h is given by: 

h e [0, 4.2] . (28) 



(b) Models where Amf^/Amf^^ ^ Arrid/ Arris: The relevant region in the p — rj plane 
is shown in Fig. ^ and the resulting allowed region in the h — a plane is presented in 
Fig. 1^. The admissible range for h is the same as in eq. (0), but the excluded area 
in the h — a plane is somewhat larger. 

(c) Models in which 5^*^ is significantly smaller than unity: The relevant region in the 
p — f] plane for this class of models is shown in Fig. |^ and the resulting allowed region 
in the h — a plane is presented in Fig. |^. The excluded region is significantly smaller 
than in scenario (a) and (b). The admissible range for h and cr are given by: 

h e [0.14, 4.18] and a E [0°, 180°] . (29) 



Unlike for scenarios (a) and (b) in this case h has also a non-trivial lower bound 
and a has an upper bound. Consequently NP contributions to the B — B mixing 
amplitude are required in this scenario. 
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IV. SPECIFIC NEW PHYSICS MODELS 



In this section we focus on two specific New Physics models that belong to the general 
class of models discussed above. We discuss first models with minimal flavor violation 
(MFV) and subsequently we present some new results relevant to LRS models with spon- 
taneous CPV. 



A. Minimal Flavor Violation 

Models with MFV comprise the SM and those extensions of the SM in which all flavor 
changing interactions are described by the CKM matrix. These models do not introduce 
any new operators implying the absence of any new CP violating phases beyond the KM 
phase. The only impact of MFV models are modifications of the Wilson coefficients of 
the SM operators due to additional contributions from diagrams involving new internal 



particles 1 10,18- 2C . 



1. Constraints from Ru, Am^^s and ek 

In the framework of this analysis the only relevant modifications of the SM predictions 
concern the mass difference Am^ of the Bd system and the parameter sk describing CP 
violation in — mixing. It turns out that for both parameters the new physics charm 
and charm-top contributions are negligible and that additional contributions to the SM 
box diagram with top quark exchanges can be described by one single parameter Fu- 
This parameter effectively replaces the relevant Inami-Lim function 5*0 (x^) appearing 
and (PD, which have to be substituted by ||18 



m eqs. 



* |2 

tb\ ' 



(30) 



and 



V 



[1 - p)r]2FuA^ + 



K 



0.226, 



(31) 



respectively. The constraint on the unitarity triangle from [c.f. eq. (§)] remains unaf- 
fected in MFV models, because it is dominated by tree-level contributions. Furthermore, 
due to the absence of any new phases in the B — B mixing amplitude the SM prediction for 
0"il)Ks [c-f- eqs. (|0])-(|^)] is unchanged. Finally, Arrid/ Arris in eq. (H) remains unaffected 
since the NP contributions cancel each other in the ratio. 

The important point we want to make is that for any given set of the parameters 



A, Bk, Pc 



and 
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one can obtain exact expressions for p, rj and Fu- Combining eqs. (|^), ( pOf ) and (|3T| ) results 
into a quartic equation that can be solved analytically. For example, for p one obtains 



{Rl-p') (l-p)CA + Ce(l-2p + C 



Ce{l 



2p + Rl] 



(33) 



where 



/ 



Glr]BmB,M^A^ 



0.226 



and 



Pc 

A^V2 



(34) 



We have used the analytic expressions to obtain the allowed regions for p, t] and Fu by 
scanning over the intervals for R^, A, B^, Pc and Aj in Tab. 2. A priori up to four solutions 
are possible, but for the present ranges of the parameters Ru, A, B^, Pc and Aj only two 
solutions are real (see Fig. ^ and Fig. |^). Adding the constraint from Am^/Ams in eq. (^ 
practically the entire solution shown in Fig. |^ is ruled out and only the second solution in 
Fig. ^ remains viable. We find that Fu is restricted to the interval 



Ftt e [1.2, 5.7] 



(35) 



Soixt) ~ 2.46 (mt/170 GeVy-^'^ 



which includes the SM value -Fj^^^ 

The value of sin 2/3 within models of MFV has obtained considerable attention re- 
cently [^,0,0 . In particular it was shown in ref . ||30| that there exists a lower bound on 
sin 2/? in these models. The existence of this bound is a non-trivial result of the correlation 
between Arriii and Ek which holds in models of MFV |3^. Using the analytic solutions for 
p, rj and Fu we have investigated the behavior of sin 2/5 as a function of Ru, A, Bk, Pc and 
Af. We find that it takes its minimal value, {sin2f3)min, when R^ is minimal and all the 
other parameters are maximal. While the size of the exact expressions we obtained for p, 
rj, Fu and sin 2/? preclude its presentation here, one can use the results to expand these 
parameters around a particular value. For example a linear expansion of sin 2/3 around 

j^min^ j^max^ Qrnax ^ prnax ^rnax gj^^g. 



sin 2/3 ~ 0.52 + 0.23 (i?„ - 0.34) - 1.09 (A - 0.86) - 0.47 (Bk - 1.00) 
-0.31 (Pc - 0.36) - 3.43 x lO'*^ ' 



^ ^ - 1.04 X loA . 

vps-i GeV-2 J 



(36) 



This expansion provides a good approximation to the exact expression for sin 2/? within 
most of the allowed parameter space for Ru, A, Bk, Pc and Af. Moreover we find that for 
reasonable changes in the boundaries of the respective intervals, sin 2/3 remains minimal 
when Ru is minimal and all the other parameters are maximal. Thus one can use eq. (^) to 
recalculate (sin2/3)mm for somewhat different values of j^i^^^^ ^ /[I'^o-^ ^ B'^'^^ , P^°-^ and A™"^. 
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For example using the more conservative values of ref. [T^ one recovers the corresponding 
lower bound on sin 2/5 which is somewhat weaker than our result, 



{sm2(3)min = 0.52 . (37) 
We note that the upper bound on sin 2/5, 



(sin2/?)„„, = 2i?rVl - (^r^)' = 0.78 , (38) 

is determined by the tangent to the Ru annulus as in the SM 

Conversely, the range of sin 2/5 due to the measured value of a^Xs restricts the parameter 
Ftt- Remarkably, taking the one-sigma range in eq. (^, only the lower bound on Fu is 
slightly improved to -F™" = 1.3. The reason is evident from Fig. Even though the 
upper bound on sin 2/5 removes a substantial part of the allowed region, the remaining 
solution (below the upper gray line) still extends over a wide range for the parameter Fu 
(as indicated by the gray-scale). 

Finally we note that the allowed interval for Fu is consistent with the result obtained in 



the generic framework discussed in section III. Indeed the ranges for h at a = or a = TC 



in Fig. ^ include the values for h'^^^^ defined via 

j^MFV^i.-'^- = (Fii/F,f ^ - 1) G [-0.48, 1.51] , (39) 

where - due to the absence of new phases in MFV models - either cr*^^^ = if Fu > Fu^ 
or a^^^^ = vr if Fu < Ff/'' . 

To conclude, models of MFV provide a consistent explanation of all the relevant exper- 
imental data. However, due to the Arrid/ Arris constraint the allowed region in the p — rj 
plane is similar to the one of the SM and Fu is close to -F^f*^. Therefore with the set of 
observables studied in this section it is very difficult, if not impossible, to distinguish MFV 
models beyond the SM. Therefore, we turn now briefly to a second set of observables which 
might help to disentangle such models from the SM. 



2. Constraints from Y Xvu 



Let us consider the decays tt^vv and B XsW within the MFV models. They 



are both related to the generalized Inami-Lim function, Xyy [jT8|,|T9[. The theoretical 
prediction for tt'^uu is given by [|Tn,|T3]: 



Bk = BR(i^+ ^ n+uu)/K+ ^ S^[A^riX,^y + [Po{X) + A\l - p)X,p]'| 



(40) 



where the values of Po{X) and are given in Tab. 1. 

Using eq. ( ^OD from the experimental value of BR(i^"'" tt'^uu) (given in Tab. 1) and 
the value of we get: 
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0.7 = S™" < {A^X.^Rt)" + 2P^{X)A^X,, (1 - p) + Po{Xf < = 11-9 , (41) 

where B"^^^ {B^°'^) denotes the minimal (maximal) value of Bk- From eq. (^) it follows 
that X^p e [X™", X^"^] or X^p e X!!^^"], where 



1 



± " Am? 



Po{X){p - 1) ± JR^BT'' - v'PoiX) 



X. 



max 



± " A^R? 



Po{X){p - 1) ± JR^BT^ - v'PoiX) 



(42) 
(43) 



A scan of the allowed p — rj region in Fig. ^ (including the a^xs constraint) and the intervals 
for A and Po{X) yields the following two ranges for (Xup): 

0.5 < X^p < 7.8 or - 9.8 < X^p < -1.6 . (44) 

Note that the uncertainty related to the measurement of BR(i^'^ tt^uu) (given in Tab. 1) 
is based on observation of a single event. Therefore, the ranges given in eq. (|4^ correspond 
to a confidence level of somewhat less than 67 %. 



Within MFV models the decay B ^ X^vv is also related to the parameter X^p [Il0| , |18 
The prediction for the branching ratio, 

BR(5 ^ X,z/z/) ^ 1.5- 10"^ -X^-, (45) 

in combination with the experimental bound (given in Tab. 1): 

BR(5 ^ Xsuv) < 6.4 • 10"^ , (46) 

provides a second constraint on X^p: 

\X,p\ < 6.6, (47) 



which excludes some of the allowed range given in eq. (^4[) . Combining the two ranges 
yields: 

0.5 < X^p < 6.6 or - 6.6 < X^p < -1.6 . (48) 

We note that in principle also the decay B Xdvv could provide a constraint on 
Xyp, but since it is suppressed by Vtd/Vts with respect to i? ^ X^z/z/ its branching ratio is 
decreased by roughly an order of magnitude and therefore at present less useful. 

Currently the constraints are too weak to rule out one of the two intervals, but the 
forth-coming data on K and B decays with final neutrinos will eventually allow to pin 
down the value of X^p and maybe provide evidence for MFV beyond the SM if X^^p turns 
out to be different from its SM value X^^ = X{xt) ~ 1.5 |T^. In this context also the 
observable {e'/e)K could be useful, since it depends on a related parameter Xuu, which in 
some models roughly agrees with Xi^p (compare for example the values for X„^ and Xj,p in 
ref. m). 
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B. LRS Model with Spontaneous CP Violation 



The spontaneously broken LRS (SBLR) model has been studied in detail in refs. [pl|-|28[] 
(see also references therein). Here we focus on the impact of the measurement of a^Ks 
eq. (H) on this model. The important feature of the SBLR model is that essentially all 
the phases and therefore the CPV observables depend on one parameter. In a certain 



phase convention [22-25] this parameter is written as rsina with |r| < m^/mj (see also 
refs. pl|j27[| ). Due to this single parameter the model is very predictive. 

The problem is that part of the model has not yet been solved analytically. In particular 
analytic expressions for the phases of the left and right CKM-like mixing matrices exist 
only within the "small phase approximation" ||2T|-p5|, which is valid for rsina ^ 0.01. 
Nevertheless a very thorough numerical analysis, beyond the small phase approximation, 
has been performed by Ball et. al. in ref. ||2^, calculating the predictions for several CPV 



observables and discussing the limitation of the small phase approximation. In particular 
it was found that within the LRSB model the measurement of Ek in eq. (|T|) and other 
observables implies 

{a^Ks)LRSB ^ 0.1 , (49) 
which is inconsistent with the combined measurement in eq. (H). It is important to note. 



however, that the analysis in ref. only used the central values for various input param- 
eters which are subject to theoretical and experimental uncertainties. Therefore we find 
it timely to reinvestigate the prediction in eq. (^) including these uncertainties and using 
the most recent values for the central values. 



1. The Analysis 



We restrict our analysis to the regime in which the small phase approximation is valid, 
which allows us to use the analytic expression for various observables of the K and B 
system. 

We scan over the entire parameter space of the relevant input parameters according to 
their uncertainties in order to find the subspace consistent with all the measured observ- 
ables. In particular, we test all the 32 different sign combinations of the quark masses, and 



allow the CKM phase 6 to be centered both around and tt (c.f. discussion in ref. [^]) 
which yields altogether 64 different signatures. For each permissible subspace we calculate 
the corresponding predicted range for a^Xs- 

Most of the analytic expression in the small phase approximation have been calculated 
in refs 



23,24 



For notations and phase conventions we refer the reader to refs 
(Note that the parametrization of CKM matrix used in ref. [E^] is the Kobayashi-Maskawa 
convention and not the Chau-Keung convention |3^, which is commonly used.) All the 



parameters relevant to our analysis can be found in Tab. 2. 
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For the kaon system we consider the observables Ami^-, Ek and {e' /e)K- The effective 
Lagrangian 

£AS=2^ that describes the AS = 2 processes, is given in eq. (19) of ref. ||23|| . The 
relevant phases of the left and right mixing matrices, 6i, 62, 7 and 6 have been calculated 
in ref. |2^. We use the following three constraints: 



(i) Am/<: Due to the lack of reliable predictions for the long distance contributions we 
require that the short distance contribution should at most saturate ArriK, i-e. 



2|Mi^| < Am^P, 



(50) 



where M^^ = {K^\C^'=^\K^) / {2mK). 
(ii) Ek- We use the upper and lower bound given in ref. P7[ : 



6.375 ■ 10^3 - 0.01 ■ 



lTeV\' 



Mo 



< 9m < 6.523 ■ 10-^ + 0.01 ■ 



ITeVy 
M2 ) 



(51) 



where, in our case. On 



2ReMi2 
AniK 



sin (argMi2) 



(iii) TZe{e' /e)K'- Due to the large hadronic uncertainties in the calculation of e'^ we only 
demand that the sign predicted by the SBLR model is positive in order to be consis- 
tent 1^ with recent measurements Thus, for the computation of 7ie{e' /e)K it 



is sufficient to use the rough approximation given in eqs. (27) and (28) of ref. 



For the B system we consider Am„ 



d,s) and a^j^g . The value of h}f^ = 



q LR 



where Mf^ = {B^\C^^=^\B'^) /{2mg), has been calculated in ref. |27 



1.7 0.051 + 0.26 In- 



Mo 



1.6 TeV 



1.6 TeVy 



Mo 



7 TeV y 
Mh ) 



(52) 



For 0"^^ an exact expression as a function of the fundamental parameters does not exist. 



so we use the small phase approximations p4 

sin^ ti 



sm a 



LR 



rjdTjb r sm a 



+ 



■ LR ■ / Vcrric , ritmt 
sm ^ risTji, r sm a \ 1 



(53) 
(54) 



where ti, t2 and denote the three angles of the CKM matrix in the Chau-Keung conven- 
tion and ?7j = ±1 parameterizes the sign associated with the quark mass mj. For our 
analysis we use the following three constraints: 
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(i) Arrid. Eq. ([T9|) , eq. (p2D and eq. (p3D yield the constraint: 



0.48 < 



< 1.05 



(55) 



where V,^ is the mixine matrix of the left handed current of the weak interaction, 
(ii) Eq. (|0D, eq. (H), eq. (H), and eq. (||) yield the constraint: 

(\/i*)2(l + /iLRe^'^rf^" 



< 0.98. 



(56) 



(iii) a^Ks- III SBLR model a^Xs is given by [ p7[] : 



a^Ks - sm 



2/5cKM + arg (l + /i^^e^' 



where 



PcKM = arg 



We require that a^Ks is positive. 



Results 



(57) 



(5^ 



We find positive values for a^Ks that are consistent with all the experimental constraints 
for several signatures. The allowed range for a^Ks that results from scanning over the 
allowed ranges of the relevant parameters in Tab. 2 is: 



a^Ks e [0) 0.3] . 



(59) 



In general, large values for a^Ks occur for large values of r sin a for which the small phase 
approximation is less reliable. Still, even for the rather large values of r sin a one still 
obtains a rough estimate for the true value of a^Xs- Therefore the fact that the range 
given in eq. ( pU]) is marginally consistent with the world average of a.^Ks sq. (^ gives 
strong indication that the inclusion of the uncertainties of the relevant parameters relaxes 
the upper bound on a^Ks i^ ^q. ( ^9]) enough to make an exclusion of the LRSB model 
premature at this stage. 
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V. CONCLUSIONS 



The CKM picture of the Standard Model (SM) is currently being tested with unprece- 
dented accuracy in various experiments. In particular the B factories BELLE and BaBar 
provide highly interesting results. There is great hope that those or future experiments 
will reveal inconsistencies that indicate contributions to flavor physics from New Physics. 

In this paper we have studied extensions of the Standard Model where the charged 
current weak interactions are governed by the CKM matrix and where all tree-level decays 
are dominated by the Standard Model contributions. We have constrained both analytically 
and numerically the ratio h and the phase difference a between the New Physics and the 
Standard Model contributions to the mixing amplitude of the neutral B system using 
the experimental results on R^, Am^ ,,, Ek and a^Ks- Fo^^ generic models we find that 
h < 4.2 and the allowed region in the h — a plane shown in Fig. |^. Models where the 
ratio Arrid/ Arris takes the Standard Model value are only slightly more constrained than 
the most general scenario of our framework (see Fig. However imposing a small CP 
violating phase significantly reduces the allowed parameter space for h {h E [0.14, 4.18]) 
and a [a E [0°, 180°]), see Fig. 0, which requires NP contributions to the B — B mixing 
amplitude. 

We have presented some new results for models with minimal flavor violation, pointing 
out that the three parameters p, rj and Fu can be obtained analytically from the experi- 
mental results on R^, Amd and Ek for any given set of rather well-determined parameters. 
Using the exact expressions we have updated the allowed interval of Fu {Fu G [1.3, 5.7]) 
as well as sin 2/3 (sin 2/3 G [0.52, 0.78]) [for which we provide an expansion in terms of 
those parameters in eq. (p6D]. Taking into account the analytical relation between p, rj and 
Fu as well as the constraints from a^Ks and Airid/ Airis improves significantly previous 



results 12 



We also consider the spontaneously broken left-right symmetric model and perform a 
numerical analysis using the "small phase approximation" in order to probe the viability 
of this model in view of the recent results for a^Xs and other observables. We find that 
the inclusion of the uncertainties of various input parameters relaxes the upper bound on 
a^Ks significantly and conclude that at present the model is still viable. 
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APPENDIX A: 



1. Global constraints on and 6^ 

To each of the three scenarios (a), (b) and (c) corresponds a different allowed area in 
the p — T) plane (shown in Fig. Fig ^ and Fig ||, respectively). Here we compute the 
global constraints on Rt and (3 for each scenario. From eqs. (|]) and (|20|) it follows that in 
scenario (a) the allowed interval for Rt{p,ri) [defined in eq. (^] is given by 

0.57 = 1 - i?;^'^^ < Rt{p, 7])<1 + i?™"^^ = 1.43 . (Al) 

In scenario (b) the upper bound on Rt is significantly stronger, namely Rt < 0.99 due to 
eq. (pOl). In scenario (c) there are two ranges, corresponding to the two allowed regions in 
the p — rj plane in Fig. ^. The upper bound on Skm implies that we have to exclude the 
interval 



0.68 = V(l - hy + hj< Rtip, v) < V(l + + = 1-33 , (A2) 

where Is = -R^*"" cos 6"^"'^ and hs = i?™" sin S"^"-^^ from the range in eq. ( [Al| ) in this scenario. 

In scenario (a) and (b) the angle /3(p, rj) [defined in eq. ([TT|)] is restricted by the maximal 
value of Ru- 

\(3{p, T]) I < arcsin(i?^"^') = 25° . (A3) 

In scenario (c) we have 

where the plus sign refers to the region with negative p and the minus sign to the region 
with positive p. Using the above equations one obtains the allowed intervals for and Od 
quoted in section |111 A| in eqs. (P2])-(B1). 



2. Analytic Boundaries in the h — a Plane 

In this appendix we describe in detail how the allowed intervals for and 6d are 
translated into the allowed region in the h — a plane. To be specific we will use the global 
extrema for and 9d in section |III A| [c.f. eqs. (p2D-(p^)]. Note, however, that essentially 



the same equations apply for the numerical analysis, where we consider each viable point 
in the p — t] plane separately. 
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a. The constraint 



In order to obtain an equation that is independent of 6d we take the absolute square of 
eq. dlBp: 



rl = ^h^ + 2/i cos a + 1 . (A5) 

Solving for h one has 



(a, Td) = — cos cr ± y cos^ cr + r J — 1 . ( A6) 

It follows that h e [/i™", /i™"^] where 

^ max(0, /.-(a, r^-)) and h^r^a) ^ h+{a, rTl (A7) 

and that h is excluded from the interval , /i^^], where 

h-^{a) ^ max(0, /i- (a, rr")) and /i+(a) = /i+(a,rr") . (A8) 

The function k^^^i^a) corresponds to the upper solid curves in Fig. |^ and Fig. ^ while the 
functions hf^{a) combine to the closed dotted curve that excludes the area around h = 1 
and 0" = vr. Due to the stronger bound on r^^*" in scenario (b) the corresponding excluded 
area is larger as can be seen in Fig. In scenario (c) there are two possible intervals for 
Trf, one corresponding to the large band between the solid curves for positive p and one 
corresponding to the area between the dotted curves for negative p (see Fig. [Zp. 

h. The 9d constraint 



To make use of the upper and lower bounds on 6d in eqs. (]2^)-(p4D we eliminate 
from the eq. (|16|), getting: 

. hsma , ^ l + Zicoscr , 

sm 29d = , and cos 29d = , ■ (A9) 



Solving for h yields 

sm a — sm Wd 

In order to find the upper and lower bound on /i as a function of o from the global extrema of 



Qd in eqs. (|2^)-(p^ we have to treat carefully the discrete ambiguities of the trigonometric 



functions in eq. ( |A10| ). To be specific in the following we discuss the constraints in scenario 



(a) and (b) that arise from eq. ( [A3D . For scenario (c) one has to consider separately the 
two allowed interval for Qd due to eq. (^^). We investigate separately the following three 
cases: 
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1. 29(1 G [~90°, 0°]: From 29^ > —32° it follows that sin29d increases monotonically in 
the interval [—0.56, 0]. Due to eq. ( |X9| ) the negative sign of sin26'rf imphes that also 
sin cr is negative and therefore a G [180°, 360°]. Taking into account that cos2^rf > 
(and therefore h < —1/ cos a if a G [180°, 270°]) we find that there is only an upper 
bound on h function of a: 

h < /i^/^(fT) = for a G [180°, 360° - |207"|] . (All) 



2. 29fi G [0°, 180°]: In this regime there are no constraints on 9^- Since sign(sin2^d) = 
sign(sin(T) it follows that also for a G [0°, 180°] we cannot derive a limit on h. 

3. 29d G [180°, 270°]: From 29d < 212° it follows that sin26'd increases monotonically in 
the interval [—0.53, 0]. Due to eq. ( |A9| ) the negative sign of sin2^rf implies that also 
sino" is negative and therefore a G [180°, 360°]. Taking into account that cos29d < 
(and therefore h > —1/ cos a if a G [270°, 360°]) we find that there is only a lower 
bound on h: 

h > Kj^[a) = h+{a, ^r"") for a G [180°, 29'^''''] . (A12) 



Together the constraints in eq. (|A11|) and eq. (|A12|) exclude the region between the 
dashed [h^^"'{(T)] and the dashed-dotted [h'^f'^{a)] curves in Fig. ^ and Fig. |^. For scenario 
(c) the analysis is similar, resulting in the dashed (p > 0) and dashed-dotted curves (p < 0) 
in Fig. 0. 
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Tab. 1: Input values 



y Q in f^T 


\I In 

V dlLlC 


ref. 


T A 

Vud 


u.y / OO 




H H II 

111 


Vub 


U.UUoOO ± U.UUUoO 




H H II 


Vcb 


U.U4U ± U.UUiO 




11 


Arrid 


yjAo ( ± U.Ui4j ps 




i 1 1I 

111 


AfriB, 


> io ps 




111 


\£k\ 


U.UUZZO 




1 1 1 

111 


D5x\K — ^ vr z/z/j 


/I I-+3.4N , , 1 n— 10 




36 




< 0.4 ■ iU 




37 


Vi 


i.do ± U.Z 




1 1 1 

111 




n PIT 




1 1 1 

111 


V3 


U.4 / 




11 


Vb 


U.OO 




1 n 1 

1C| 


Vx 


n QQ/1 




ml 


s 


1 14 ± 05 




111 


Pc 


0.31 ±0.05 




m 


PoiX) 


0.42 ±0.06 




m 


rfitimt) 


(166 ± 5) GeV 




m 


A 


0.819 ±0.040 




111 


A 


0.2237 ±0.0033 




111 


Bk 


0.87 ±0.13 GeV 




111 




0.23 ± 0.025 GeV 




111 




4.11 X 10-" 




TOfl 
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Tab. 2: Input values for the SBLR model 





\/ 1 n #3 
V dl Lie 


1 Ci. 


mu[Mz) 


U.UUzo ± U.UUU4/ (oreV 






U.UU4/ ± U.UUUo LreV 




mc[Mz) 


U.D^ ± U.UOD VoreV 


1 9C 1 

|38| 


ms{Mz) 


U.Uyo4 ± U.UiZ voreV 


1 9C 1 

|38| 


mt[Mz) 


1C1 _i_iQ r^^A/" 


1 9C 1 

|38| 


mb[Mz) 


6 ± U.ii VjreV 


1 9C 1 

|38| 




snnn + 5nnn GpV 




Mh 


12000 ± 2000 GeV 


|27| 


Sk 


0.35 ±0.1 


|22| 




1.2 ±0.2 


|27| 


ti 


0.231 




t2 


0.041 ±0.018 




t3 


0.016 ±0.0016 




V4 


5± 1.5 


P 
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FIGURES 
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FIG. 1. Determination of the apex {p,rj) of the unitarity triangle: The grey area is the 
presently allowed region, which is determined by (a) the measurement of ll^fe/V'cbl (dashed cir- 
cles), (b) the observed — mixing parameter Am^ (dotted circles), (c) the lower bound on 
the B^ — Bg mixing parameter Am^ (dashed-dotted circle), (d) the measurement of ex (solid 
hyperbolae) and (e) the combined result of the CDF, BELLE and BaBar measurements of a^Xs 
(thick grey lines). 
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V 




FIG. 2. The shaded area corresponds to the allowed region in the p — rj plane for scenario (a). 
It is constrained solely by i?„ G [0.34, 0.43], which is determined by the measurement of |T4fe/V"c6| 
(dashed circles). For comparison we also show the other constraints appearing in Fig. ||, which 
are not relevant here. 
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FIG. 3. The dark area corresponds to the allowed region in the h — a plane for scenario 
(a) according to our numerical analysis. It is embedded in the light area that results from 
the analytical boundaries based on the global constraints on r^i from Am^^g and 6d from a^Xs- 
Specifically the allowed range of excludes h to be above h'^^^{a) (solid curve) or between hf^{a) 
(dotted curves), while the allowed range of 0^ exclude the region between h^^^{a) (dashed curve) 
and /i^"^ (cr) (dashed-dotted curve) . The numerical results are not much more constraining than 
the analytical ones, except for the "holes" in the dark regions that correspond to the excluded 
region within the Ru annulus, which has been ignored in "global" analysis. (See appendix |A 2| 
for details.) 
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FIG. 4. The shaded area corresponds to the allowed region in the p — r] plane for scenario (b). 
It is constrained by Ru G [0.34, 0.43] (dashed circles) and Ams/Anib (dashed-dotted circle). For 
comparison we also show the other constraints appearing in Fig. |l], which are not relevant here. 
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FIG. 5. The dark area corresponds to the allowed region in the h — a plane for scenario 
(b) according to our numerical analysis. It is embedded in the light area that results from 
the analytical boundaries based on the global constraints on from Amj^ and from a^^s- 
Specifically the allowed range of excludes h to be above h^^^{a) (solid curve) or between 
hf^{(j) (dotted curves), while the allowed range of O^i excludes the region between h^^'^{a) (dashed 
curve) and h'^^^^a) (dashed-dotted curve). Due to the stronger bound on rj**" in scenario (b) 
the excluded area from hf^{a) is larger than in Fig. |^. (See appendix A 2 for details.) 
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FIG. 6. The shaded area corresponds to the allowed region in the p — r] plane for scenario (c). 
It is constrained by Ru G [0.34, 0.43] (dashed circles) and tan 6km < 0.25 ("small phase"). For 
comparison we also show the other constraints appearing in Fig. |l], which are not relevant here. 
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FIG. 7. The dark area corresponds to the allowed region in the h — a plane for scenario 
(c) according to our numerical analysis. It is embedded in the light area that results from 
the analytical boundaries based on the global constraints on from Arud^s and 6^ from a^Xs- 
Specifically there are two possible intervals for r^, one corresponding to the large band between 
the solid curves for positive p and one corresponding to the area between the dotted curves for 
negative p. The allowed ranges of 0^ limit the allowed region to be between the dashed {p > 0) 
and dashed-dotted curves {p < 0). Note that unlike for scenarios (a) and (b) in scenario (c) h 
has a lower bound and a has an upper bound. (See appendix |A 2| for details.) 
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FIG. 8. The shaded area corresponds to the "first" solution resulting from eqs. (^), ( |30[) 
and (^) within models of MFV. The grey scale indicates the magnitude of Fu as a function 
of p and r]-. The lightest region corresponds to Fu — 0.4 and the darkest region corresponds to 
Fft — 2.7. Adding the constraint from Am^/Amb (dashed-dotted circle) practically the entire 
solution presented in this figure is excluded. 




FIG. 9. The shaded area corresponds to the allowed region according to the "second" solu- 
tion resulting from eqs. (^), ( pOD and (^) within models of MFV. The grey scale indicates the 
magnitude of Fu as a function of p and r]: The lightest region corresponds to Fu — 1.2 and the 
darkest region corresponds to Fu — 5.7. 
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